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1 Introduction 



One of the most studied sectors of A/" = 4 SYM is probably the so-called s/(2) sector. In the planar 
limit, i.e. number of colours N ^ oo and SYM coupling qym — so that the 't Hooft coupling 

A = 2^?^j,,iV = Stt V (1-1) 

stays finite, it results from local composite operators with the form 

Tr(P"Z^) + ...., (1.2) 

where V is the (symmetrised, traceless) covariant derivative acting in all possible ways on the L 
bosonic fields Z. The spin of these operators is s and L is the so-called 'twist'. Moreover, this 
sector would be described - thanks to the AdS/CFT correspondence |1] - by string states on the 
AdSs X spacetime with AdSs and charges s and L, respectively. In addition, as far as the 
one loop is concerned, the Bethe Ansatz problem is equivalent to that of twist operators in QCD 
[21 [3]; and this partially justifies the great interest in this sector (11. 2p . 

In general, the eigenvectors of the dilatation operator are suitable superpositions of the opera- 
tors (II. 2p with definite anomalous dimension (eigenvalue) 

A{g,s,L)=L + s + ^{g,s,L), (1.3) 

where ''y{g,s,L) is the (renormalised) anomalous part. According to the AdS/CFT strong/weak 
duality, the set of anomalous dimensions of composite operators in A/" = 4 SYM coincides with 
the spectrum of energies of the corresponding string states, although the perturbative regimes are 
interchanged. The highly nontrivial problem of evaluating anomalous dimensions in AA = 4 SYM 
was greatly simplified after the discovery of integrability in the purely bosonic so(6) sector at one 
loop [1]. Later on, indeed, this fact has been extended to the whole theory and at all loops in the 
sense that, for instance, any operator of the form (11.21) is associated to one solution of some Bethe 
Ansatz-like equations and then any anomalous dimension is expressed in terms of one solution [5] . 
Actually, this is only part of the full story, although the rest should not worry us in this context. 
In fact, an an important limitation emerges as a consequence of the on- shell character of the 
asymptotic Bethe Ansatz: as soon as the interaction reaches a range greater than the chain length, 
then it becomes modified by unpredicted wrapping effects. More precisely, the anomalous dimension 
is in general correct up to the L — 1 loop in the (SYM) convergent, perturbative expansion, i.e. 
up to the order g'^^~'^. Which in particular implies that the asymptotic Bethe Ansatz should give 
the right result whenever the subsequent limit (II. 4p is applied, as in the leading expansion (11.50 . 
The important large twist and high spin limit goes as a double scaling as follows: 

s— *oo, L^oo, i = - — = fixed . (1.4) 

m s 

Yet, the string theory becomes dual as semiclassical expansion with the string tension — +oo 
playing the role of inverse Planck's constant. This means that this limit is to be considered before 



2 



the scaling (11.41) (cf. for instance ^ and references therein), thus imposing a different hmit order 
with respect to our gauge theory approach (cf. below for more details). The relevance of this 
double limit has been suggested in first instance by [7J within the one-loop SYM theory and then 
motivated in [6] and in p| within the string theory dual. In fact, calculations of the latter authors 
pointed towards the following generalisation (at all loops) for the the anomalous dimension formula 
found in [Zjl] 

7(^?, s, L) = /((?, j) Ins + 0((lns)-°°) . (1.5) 

Moreover, by describing the Bethe Ansatz energy through a non-linear integral equation (like in 
other integrable theories [10]), this Sudakov scaling has been remarkably confirmed in p]. There 
this statement was argued by computing iteratively the solution of some integral equations and 
then, thereof, the generalised scaling function, f{g,j) at the first orders in j and g^: more precisely 
the first orders in g"^ have been computed for the first generalised scaling functions fn{g), forming 
the crucial expansion (see below for the motivation) 

oo 

f{9,j) = J2f-(9)r- (1.6) 

n=0 

As a by-product, the reasonable conjecture has been put forward that the two-variable function 
f{g,j) should be bi-analytic (in g for fixed j and in j for fixed g). Concentrating on the fn{g) one 
easily realizes that in order to check the AdS-CFT correspondence the knowledge of their behaviour 
for all values of g is extremely important. In particular, comparisons with string theory results 
involve the strong couphng limit of all the /n(fi')'s, after the papers [IIl|12l[T3] clarified the large g 
asymptotic expansion of fo{g) = f{g)- In this context, in paper [H] it was shown how to get the 
contributions beyond the leading scaling function f{g), by means of one linear integral equation, 
which does not differ from the so-called BES equation, which covers the case j = 0, cf. the last 
of [5], but for the inhomogeneous term. Re- writing this equation in the form of a linear system, 
in the paper [15| the computation of the first generalised scaling function fi{g) at strong coupling 
was performed, showing its proportionality at the leading order to the mass gap m{g) (see f l3.13p of 
the 0(6) nonlinear sigma model (NLSM) and giving, consequently, the first evidence of the related 
proposal^ of Alday and Maldacena [H|. Afterwards, the embedding of the 0(6) NLSM into Af = A 
SYM at large g was shown on formal grounds in [16], where also the leading contribution in the 
strong coupling limit of fs^g) was computed. In our contemporaneous paper |17], we could write 
a set of linear systems of integral equations (of the Fredholm type), whose solutions could give 
all the generalised scaling functions. The n-th system, which gives the n-th generalised scaling 
function, needs information involving the solution of the m-th system, with m < n — 3 and this 
from the practical point of view allows us an iterative solution. Interestingly, important structures 
describing the fnid) as functions of g started to be identified as the values of the density of Bethe 

^0{\ns~°°) means a remainder which goes faster that any power of Ins: hm (In s)'"'0((ln s)^°°) = 0,Vfc > 0. 

s — >oo 

^According to this proposal, formulated on the string side of the correspondence, when g oo and j /m{g) is 
fixed, f{g,j) + j should coincide with the energy density of the 0(6) nonlinear sigma model. Results in the present 
paper and of [IHIIII] concern the case j <C m{g). 
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roots and holes and its derivatives at the origin. Exphcit evaluations at large g could be done 
and the strong coupling limit of fsig) and f4{g) was explicitly worked out, finding agreement with 
analogous results [18] from the 0(6) nonlinear sigma model. As stressed in [17], we remark that 
the agreement for f^lg) is highly nontrivial, since it involves the full interacting theory of the 0(6) 
NLSM. 

In the present paper we want to present a systematic approach which allows us to tackle the 
computation of all the generalised scaling functions /„((?) and consequently of f{g,j). 

The plan of the paper is as follows. In section 2 we will write general equations for the densities 
of roots and for the anomalous dimension, which are valid in the regime of large spin s, upon 
neglecting terms 0((lns)~°°). In the limit (11. 4p such equations will be rewritten (section 3) as an 
infinite series of linear systems of the Fredholm type, whose solutions give the complete (i.e. all j) 
generalised scaling function. In Section 4 we will develop a systematic technique which allows to 
write the solutions of the n-th system - which gives fn{g) - in terms of the solutions of 'reduced' 
systems and of the values of the density of Bethe roots and holes and its derivatives in zero. In 
Section 5 the large coupling limit of all the relevant equations is studied and an explicit, recursive 
way of finding their solution is given. As an application, in Section 6 the explicit computation of 
f3{g), . . . , fsig) at strong coupling is performed. 

2 High spin equations 

In the framework of integrability in A/" = 4 SYM, it was useful [19| to rewrite the Bethe equations 
as non-linear integral equations [IDJ. In particular, this approach was pursued for the sl{2) sector 
of the theory. In this sector states are described by s Bethe roots, concentrating on an interval 
[—b, b] of the real line, and L 'holes'. For any state, two holes lie outside the interval [—6, b] and the 
remaining L — 2 holes lie inside this interval. For the lowest anomalous dimension state (ground 
state) the L — 2 internal holes concentrate in the interval [— c, c], c < b, and in this interval no roots 
are presenile]. The non-linear integral equation for states of the sl{2) sector involves two functions 
F{u) and G{u,v) satisfying linear equations. It is convenient to split F{u) into its one-loop Fq{u) 
and higher than one loop F^{u) contributions and to define the quantities <Jh{u) = -^F^(u) and 
cro{u) = -^Fq{u), which, as the spin gets large, acquire the meaning of, respectively, higher than 
one loop and one loop density of both Bethe roots and holes. 

In this context, in the large spin limit the ground state anomalous dimension depends on the 

■^c is technically known as 'separator'. 
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densities as follow^, 



bo 



dv 



-bo 



X^\V] X {V) 

i 



+00 



dv 



X'^[Vj X [V 

i 



[ao{v) + crniv)] 



:2.ii 



x^(v) X (v) 



aH{v) + 0{{\nsr'^), 



where we introduced the function Xc{u) which equals 1 if — c < u < c, where the internal holes 
concentrate, and otherwise. This means that, as far as the computation of the generalised scaling 
functions is concerned, one can rely on (12. ip . The terms depending on ao{u) get more manageable 
after using the important relation introduced in [1] 



bo 



dvf{v)ao{v) 



bo 



dvf{v)a'Q{v) + Oi(\ns) 



(2.2) 



where the Fourier transform of the function crQ(f ) satisfies the integral equation. 



a^.ik) = -4vr- 



^_e-^cos^ e-^ 



V2 



2 sinh ^ 2 sinh 



rfMe*^"Xco(M)c^o(«) - 4vr5(A;) ln2 , (2.3) 



with the parameter cq such that the normalization conditioiifl 

duxMa-oiu) = -27r(L - 2) + 0((lns)-°°) 



(2.4) 



holds. It is convenient to rewrite the previous relations entirely in terms of Fourier transforms: 



-An 



f ^ cos jl 
2 sinh 



e-V .+00 



sinh 



|fc| 



,sm{k - h)co ^ „ 



dk ^ .sin kco 



2ti 



k 



-27r(L-2) + 0((lns)-°°) 



(2.5) 
(2.6) 



These relations have to be solved together and give, for any values of L, cq and o"q(/c) at large spin 
s. In particular, we can go to the limit (11. 4p . In this case crQ{u) and cq expand as. 



Co 



[2.7) 



n=0 



n=l 



^The parameter &o > c is the one loop contribution to b: therefore, it depends on s through the solution of the 
linear equation for Fq{u) (see [U). 

^The physical meaning of (|2.4p is that, for — co < u < cq, (Tq{u) approximates the density of holes, that in the 
one loop theory fill the interval [— co, co], where no roots are present. 
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and it is not difficult to give, for instance, the values of the ffist two coefficients of the expansion 
of c: 



.(1) 



TT 

4' 



-^ln2. 

4 



On the other hand from (4.10) of [H] we get that in the large spin limit the higher than one 
loop density satisfies the linear integral equatioii^. 



1 + 



aH[u) 



—iL— In 
du 

d 



2x~{uY 



1 + 



2x+{uY 
1 



+ 00 



dvxc{v) 



-In 

du 



2x+(u)x (v) 



1 - 



+ 



2x (u)x'^(v) , 



+ i—e{u,v)+i— 

du 1 + [u — v)^ 



d 

dv— 

oo 



du 



In 



dv 



1 



2x+ {u)x~ (v) 



2x {u)x+{v) 

'•+°° dv 



TT 1 + (u — vy 

1 



+ ieiu,v) a'o{v) + 

TT 1 + [U — VY 



(2.9) 



^o(^) - 



d 

dv 

oo 



du 



In 



2x+{u)x-{v) 



1 



+ iOiu, v) 



2x {u)x'^{v) 

which has to be solved together with the conditions, 

»+oo 



aH{v)+0{{\ns) 



/ + 00 
duXciuMi^) + aniu)] = -27r(L - 2) + 0((lns) 
-oo 



(2.10) 



+ 00 



duxMa'oi^) = -2n{L - 2) + 0((lns) 



As in the one loop case, it is convenient to rewrite, in terms of Fourier transforms, equation (12. 9p . 



sinh y- 



+ 



sinh ^ 



dh 
]h\ 



r=l 



+ sgn(/i) J2 Yl Cr,r+i+2.(-l)"+"e-^ (^Jr-l{V29k)Jr+2.{V2gh) 



e 2 + 



oo 
oo oo 



(2.11) 



r=2 u=0 
Jr-liV2gh)Jr+2u 



_1M 
e 2 



sinh 



|fc| 



dp . . / sin(/c-p)c e 2 



dp . . / , , ^sin(/i — p)c 
2tt h — p . 

'~ "^^ dp ^ sin(A;-p)co 



+ 0((lns) 



^The quantity 6{u,v) appearing in (|2.9p is the well known dressing factor. For related notations, we refer to the 
last of [5]. 
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(2.12) 



and also the normalization conditions, 

^m^^^^^)^ = -2vr(L-2) + 0((ln.)-), 

-[a^,{k) + an{k)f-^ = -27r(L-2) + 0((lns)-). 
Again, in the limit (11.41) both c and (Jh{u) expand in powers of j, 

oo oo 

aniu) = [Y,^H\^)r] + 0((ln.) — ) , c = J^c^-^j^ , (2.13) 

n=0 n=l 

SO that we find convenient to use the all loops densit}0 

aiu) = aniu) + a-^iu) , (2.14) 

which, in the limit (11.41) expands as 

oo 

a{u) = [^a(")(w)j11ns + 0((lns)-°"). (2.15) 



n=0 



As in the one loop case, it is easy to give the first orders in the expansion of c in the limit (II. 4p : 

c<"^^5^. c'-'^-. ^'"^-f(°> . (2.16) 

4-ag)(0) [4-af^^n\^2 



'H 

Finally, by comparing (12. ip with (12. lip , we realize that, very easily, when s — >■ oo, 

7((7, s,L) = - hm anik) + 0((ln s)-°°) . (2.17) 

And this relation extends the validity of the analogue in [20], valid only at order Ins, to all orders 
(lns)~™, m G N. Equality (I2.17P implies, very simply, that 

f„ig) = UPiO). (2.18) 

A systematic approach to the computation of fn{g) using the solution of (12.111) and explicit formula 
which allow their exact determination at strong coupling is the main issue and the main result of 
this paper. 



^The one loop quantity crQ(u) is an approximation, according to (j2.2p . to the real one loop density. Since in this 
article we will always neglect terms 0((lns)~°°), such approximation is completely justified. 
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3 On the calculation of the generalized scaling functions 



From result (12.181) we realize that the generalised scaling functions fn{g) can be extracted from 
the n-th component a^^\u) of the solution of (12.91) in the limit (11.41) . We are therefore going to 
analyze in a systematic way (12. 9p or (12. lip in such a limit. Equations for cr^H{u) and (J^^\u) were 
already studied: the former is the well known BES equation (last of [5]), the latter was treated in 
detail in [15]. Since in this paper we will use results involving (t^hI''^) ^^"^ ^h\'^)^ "^^ "^ill treat 
also briefly these two cases, which, in addition, need to be considered separately from the rest of 
the a^^\u). 



3.1 The BES equation 

If we restrict (12. lip to the component proportional to In s we get the BES equation for a^j^{k). We 
now briefly describe how to rewrite such equation in a form suitable for our future manipulations. 
We first define 

and then expand S'^'^\k), > 0, in series of Bessel functions, 

S(»>(i) = f sg'(9):^Si^ + f sgL.to):^^^!^. (3.2) 

p=l p=l 

On the coefficients S^\g) the BES equation implies the linear system, 

oo oo 
m=l m=l 

(3.3) 

oo oo 

S^J_M = 2V2gS,,r - 2(2p - 1) J] Z2p-iM9)Si%) - 2{2p - 1) J] Z2p_i,2™-i(^)41i(^) , 

m=l m=l 

where we introduced the notation 

, f+'^dh UV2gh)JUV2gh) 
Zn,M - ^ ^^3^ . (3.4) 

The cusp anomalous dimension can be extracted from the relations 

\\m S^'\k) = ^ig) , fig) = V2gS{'\g) , (3.5) 

and its strong coupling behaviour was completely disentangled in [12]. 
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3.2 On the first generalized scaling function 



The first generalised scaling function was studied in [15]. Here we briefly recall the main results. 
We define the even function 

2 sinh ^ 



^^^nA:) = ^^^^2^(A:), (3.6) 



and introduce the two functions 



— Jr{V2gh) ar{g)= -—Jr{V2gh) ^. (3.7) 



Expanding, for /c > 0, in a series involving Bessel functions, 



p=l p=l 



(3.8) 



the coefficients Sr^^ {g) satisfy the hnear system, 

/ oo oo 

Sil\g) = 2 + 2p^-a,,{g)-2J2Z2pM9)S^l{g) + 2Y,Z,,,2n^-l{g)S^J,_,{g) 



m=l m=l 



^(1) ^^-j oo oo 

2^"^ ^ = -a2p-l(^) -2^Z2p_i,2m(^)42(^) -2 5]]^2p-l,2m-l(^)52m_l(^) 



(3.9) 



2p _ 

m=l m=l 

In paper [15] we found the following asymptotic strong coupling solution to the system (13. 9p : 



■^^-^ 1 (m — + 1) (7^" ^ 



(3.10) 



Si%) = -2mV(-l)"'-Ql^^±^%; 
2m vy; ; r m - n' + 1 o^n' 

n'=i ^ I y y 



where the coefficients 

n' 

b2n' = 2-"'(-l)"'E 



7-1 r)2A; 



{2ky.{2n' - 2ky. ' 



(3.11) 



"'-1 p o2fc 



62n'-l = 2-"'+i(-l)"'-l5^ 



(2fc)!(2n'-2A:-l) 



with i?2fc the Euler's numbers, sum up to the generating function 

oo ^ 

b{t) = V = ^ . (3.12) 
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In addition, the behaviour 

h{g) = 2 \im S^'\k) = V2gS['\g) = -l+m{g), m{g) = -^g'^e-^^' [1 + 0(1/ g)] , (3.13) 

where m{g) agrees at the leading order with the mass gap of the 0(6) nonhnear sigma model, 
expressed in terms of parameters of the underlying AdSs x sigma model, was shown after 
numerically solving the system ( 13. 9p . 



3.3 On the second and higher generalized scaUng functions 



We now give a general scheme for tackling the problem of computing the n-th generalised scaling 
function fn{g) for n > 2 at arbitrary value of the coupling constant. 
We start from (12.111) and define the function S{k): 



2sinhi^ p-^ /'+°° 

Ins S{k) = —-y^ [cTH{k) + ao^(fc)] + 



2TT\k\ 



7c\k\ 



dp r-'^^ \ - ^ ,,sm(k-p)c 1 
^Kip) + Mp)] , (3.14) 



which, differently from the cases n = 0, n = 1, depends on the all loops density (12.141) a{k) 
anik) + (To(A;). We go to the limit (fL4|) . We have 



S{k) = s^'>\k) + s^'\k)j + J2S^''\k)r- 



(3.15) 



n=2 



and we concentrate on S^^\k), with n >2. For such functions the following equation hold: 



S^''\k) 



TV 



^^^y{-^r'JAV2gk)UV2gh) 



1 — sgn{kh) 



e 2 + 



oo oo 



+ sgn(/i) J2 Yl Cr,r+i+2.(-l)''+'^e-^ (^Jr-i{V2gk)Jr+2.{V2gh) 



(3.16) 



r=2 u=0 

Jr-i{V2gh)Jr+2u{V2gk) 



_K\h\ 

sinh 



\h\ 



e 2 



sinh 



dp ^ sm(h — p)c 



271 



h — p 



where the simbol \jn means that we have to extract only the coefficient of j" in the limit (II. 4p . 
after having removed the overall factor Ins. 

Again, if we restrict the domain to > we can expand in series of Bessel functions. 



p=i p=i 

in such a way that the n-th generalised scaling function is expressed as (I2.17|3.14p : 

Ug) = V2gS{''\g) . 



(3.17) 



(3.18) 
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After some computations, we find the following system of equations for the coefficients of S^'^\k), 
with n > 2, 



m=l m=l 

oo 

Si;U9) = 4?-i(^?)-2(2p-l)5^Z2,-i,2™(^)5tH^)- (3.19) 

m=l 

oo 

- 2(2p-l)^Z2,_i,2^_i(^)^(li((7), 

m=l 

where the 'forcing terms' Ar"\g) are given by: 



.(„), , /"^°° dh Jr(V^gh) f^°^ dp sm(h — p)c^^ , ^ , 



(3.20) 



These systems have all the same kernel, which coincides with the BES one, and differ only for 
their forcing terms. The inforcing of the normalization conditions in fl3.20p will show how the n-th 
forcing term depend on the solutions of the m-th system, with m < n — 3, allowing, therefore, their 
iterative solution. This will be the subject of next section, where we are going to systematically 
tackle the problem of finding Ai"\g) for all values of n, up to the desired order. 



4 Systematorics 

The main obstacle to obtain a fully explicit expression for the infinite linear system at a generic 
order n is the double expansion contained in the term sm{{h — p)c{j)) of equation fl3.20p . A similar 
structure is also present in the normalization conditions (I2.12p . 

In order to overcome this technical problem, it is worth to remember a standard result of 
combinatorics known as the Faa di Bruno's formula [21]. Let f{x) and g{x) be a pair of functions 
admitting (at least formally) a power expansion of this kind 

oo „ oo 
n=l n=l 

then the composition g{f{x)) admits the following power expansion 

°° I, 

^7(/(x)) = Mx) = 5^^x^ (4.2) 

n=l 

where the coefficients /i„ have the following form 

n 

hn = '^gkBn,k{fl, ■ ■ ■ , fn-k+l) ■ (4.3) 
k=l 
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Bn,k{f) is the Bell polynomial defined as 

{ji,...J„_fc+i} m=l 

and the sum runs over all the non negative j's satisfying the conditions 

n—k+l n—k+l 

^ 3m = k, ^ mjm = n. 

m=l m=l 

The previous equation will be our main tool in the remaining part of the section. It is straightfor- 
ward to apply the previous formula to the present case sin(pc(j)) with 



oo 

sm X 



n=l 



en = ^r+^((-ir-l) (4.5) 

and 

(j) = E^^"^^'"- (4.6) 



oo 

c{j) = ■ 

n=l 

We end up with (we divide by p for future convenience) 

^mMj)) ^ vS^5„_,(pe«,...,Hn-A: + l)!c("-'=+i))j- = 
p ^ — ^ pnl 

n=l 

oo 

n=l 

Let us now use this result in order to write in a more convenient way both the normalization 
condition fl2.12p and the forcing term (13.201) . We begin with the analysis of A„(p) in order to put 
it in a more suitable form. Some elementary manipulations bring it to the form 



A„(p)=5^a/?n,.(c(^\...,c("-'=+^);p), (4.^ 
k=i 

with 



/5.,.(c«,...,c("-'=+^);p)= Yl (P"~') n 



{jl,...Jn-fe + l} "^=1 



n—k+l 71— fc+1 



Y jm = k, ^ mj^ =n. (4.9) 

m,=l m=l 

The forcing term and the normalization conditions have a common structure 

^'"<'" ^<^»-»(p.) = f f; AM) A ( t *'"'(P.) J-") In . = 



Pi 



^n=l / \n=0 

oo 



J2^n{pi,P2)rins, (4.10) 



n=l 
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where 

n 

r.(pi,P2) = 5^Afc(pi)a("-'=)(p2). (4.11) 

fc=l 

We can finally write the coefficient in the expansion in powers of j of the integral over the momen- 
tum which appears in the forcing term as 

/r„(/i-p,p)lns, (4.12) 

along with the normalization condition 

dp 



T^{p,p) = -7r5n,i. (4.13) 

Our next step will be to enforce the normalization condition in the forcing term in order to 
gain a simplification of its structure. We notice that A„(p) has a momentum independent term 
corresponding to the term = 1 in the sum ( 14. Sp and hence rn{pi,P2) admits the following 
decomposition 

n 

rn{pi,P2) = rf (P2) +f„(pi,p2), ri°)(p2) = J]a("-'=)(p2)c(^^ (4.14) 

k=l 

The normalization condition then becomes 

rW(p) = ^ T^ip,p) + vr5„,, , (4.15) 



) 27r 2tt 
which allows to subtract the ri^"* (p) contribution in the forcing term: 

/■+°° dv 
2j ^T„,{h-p,p) = 

r^"^ dv 

= -2n6n,i + 2 j £ [r„(/i-p,p)-r„(p,p)] (4.16) 
= -27i6n,i + 2 ^ [Anih~p,p)] . 

(2) 

A brief inspection allows to realize that such a subtraction is responsible for the fact that Ar = 
which imply f2{g) = 0, as already noticed in the literature. 

The final step is to make such a subtraction explicit for any n. In order to achieve this, we 
have to take into account the properties of parity of the Fourier transform of the densities of Bethe 
roots (they are even with respect to the momentum p), and the following relation which allows to 
perform the integral over p, 



du^ 



dp s ^{n-k)f.^\ J 1 



/a^"-^-^(p), 4 = ^(l + (-l)^), (4.17) 

u=Q ./_oo 2 



13 



which is different from zero only for s even, due to the parity properties of a*-" (p) . 
It is then possible to rewrite A„(/i — p,p) as follows 

n 

A^{h-p,p) = Y,^'^''-'\p)iMh-p) - Mp)) , (4.18) 



k=l 



where 



k k-l + l ( Jm)\jm 

Ak{h-p)-Ak{p) = J2^i Y: n ^-rr i(h - pY-' - p'-') 

«=1 {jl,-,jk-l+i} m=l 



'=1 {ji,-dk-l+l} m=l \s=0 ^ / 

^tc. E 'ff^(E(':n'''-'-(-i)v.(-P) 

m=l m=l 

The last step comes from the fact that it is always (— — p'^^ = 0, because is non-vanishing 
only for odd I. One can also notice that the subtraction and the fact that ^2 = allow the sum 
over / to begin from Z = 3. 

The previous result, together with eq. (14.171) . allows to write down, for n > 2, 

/^°°^r„(/.-p,p) = 2^''°°^ [T^{h-p,p)-T^{p,p)] (4.20) 
dp 



' — oo 

r>+00 



[An{h-p,p)] 

oo 

n k / 1-2 y , ^ ^ \ k-l+l 



k=\ 1=3 \s=0 V / / {ju...,Jk-i+i} ™=1 

k-l+l k-l+l 

m=l m=l 

which is nothing but the explicit n-th term of the j expansion of the integral over p which appears 
in the forcing term. Then, if we pose 

ji,s_^ dh Jr{V2gh) 

' Jo 27rh sinhf " ' 
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we can explicitly write down the generic expression for the forcing term Ai"'\ n > 2, entering the 
system f l3.19l) for the n-th term of the j expansion of the function S{k): 

Ai-\g) = (4.21) 



k=l 1=3 \s=0 ^ / / {h,...Jk-l+i} m=l 

k-l+1 k-l+1 



m=l m=l 



Because of the particular form of the forcing terms, it is convenient to write the solution of (13.191) 



as 



Sl!^\g) = (4.22) 

-^i:i:«.fi:(':n4(-o-v'»->^(^'fe)) e n""""*'^ 

fc=l 1=3 \s=0 ^ ^ / s-- - 1^ "i^i 

-l+l k-l+1 

3m = 1, ^ mjm = k 



fc=l 1=3 \.=0 ^ - ^ / m=l 

k-l+l k-l+1 



m=l m=l 

where the "reduced" coefficients Si'''^ satisfy the equations 



m=l m=l 

oo 

SiUg) = 4?-i(^7)-2(2p-l)5^Z2,_i,2™(^7)^S(^)- (4.23) 

m=l 

oo 

- 2(2p - 1) 5^ Z,,^,,,^^,{g)sS,_,{g) , 

m=l 

with the reduced forcing terms, 



-oo 



dh,2k-iJr{V^gh) 



which are 'known' functions, i.e. they do not depend on the quantities cr'-^'-''''^^ We notice that 
inside the structure of the forcing term aI^-' (g) (I4.2ip we find the constants c'-™'-' , with m < n — 2 
and the densities of the Bethe roots at m = (together with their derivatives) cr^"'')'^^), with 
n' < n — 3. In addition, the constants c*^™-* can be related to cr*^"'-''*^*-', with ra' < m — 1, by means of 
the normalization condition f l4.13p . thus leaving the forcing term Ai'^\g) as dependent on cr("')'(''), 
with n' < n — 3. Let us now show this. We first of all notice that, for n = 1, we have 
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and that, for m > 1, the normahzation condition takes the form 



m k k-l+l , (m')\j i 

E E (.)-'+V(-^')^('-) E n = ' (4-26) 

k=l 1=1 {ji,-,jk-l+l} m'=l 

k-l+1 k-l+1 

m'=l m'=l 

After a brief inspection of the latter it is possible to realize that, at order m, the only term which 
contains c*^'"^ can be singled out by taking k = m, I = 1, and that the remaining terms in the sums 
only contain c'-'^^ with k < m. 

As a consequence, we can write a recurrence relation 

"^-1 (m-fc);(0) 



(m) ^ _ >p ^ (fc) _ 
(^(0);(0) 

fc=l 



^{m-k);{l-l) fc-'+l 

-EE&W-^^STo^ E rn>l, (4.27) 

m'=l m'=l 

which, together with the initial condition (14.251) . allows to express all the c*^™^^ recursively, in terms 
of cr''"'-''*-'^^ with n' < m — 1. Therefore, we conclude that the forcing term Ai"\g) (14.211) and the 
solution Si^^g) (14.221) actually depend only on cr("'')'('^), with n' < n — 3, i.e. on the solutions of 
previous systems. Consequently, at least in principle, the solution for the Sr"'^ may be written in 
a recursive form. 

To summarize, the principal result of this section is formula (I4.22p : the evaluation of the n-th 
generalised scaling function /„((?) = V^gS^i^\g), for n > 2, is eventually reduced to the knowledge 
of S[''\g) and of the densities and their derivatives in zero, cr("')'(*) (|4.17p . with n' < n — 3. In 
next subsection, we will show that S[''\g) (and fi{g) = V^gS^j^\g)) can be given an integral 
representation in terms of the solution of the BES equation. However, this connection to the BES 
equation (true, for obvious reasons, also for cr(o)'('^)) is not true for the densities and their derivatives 
at zero a^'^'^'^^\ n' > 1: in order to find them, one needs more additional information, i.e. the full 
solution Sr^\g), Sr''\g), for all r, to the systems (13. 9p . (I4.23p . However, we again stress that, due 
to iterative structure of (I4.22p . an explicit solution for the Si^\g) is can be found by recursive 
methods. This will be explicitly shown in the strong coupling limit (Section 5). 

4.1 Mapping the reduced systems to the BES equation 

As stated before, the main point of this subsection is to write down an integral representation for 
the reduced coefficient s[''\g) and for S^\g), in terms of the solution of the BES equation. 
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As a first step we rewrite the BES linear system (13.31) introducing the even/odd Neumann 
expansioij^ 



a^^\V2gt) = Si'jig) J2A^9t), a^'\V2gt) = ^ Sf^M J2,-i{V2gt) , (4.28) 
p=i p=i 



with the coefficients Sr^\g) given by 



4?(^)=2(2p) 



jaf{t)J,,{t), ^£l,(^)=2(2p-l) 



dt 



a^'\t)J2p-i{t) . (4.29) 



Then, the BES hnear system can be cast in the form 




J2p{V2gt) = Q^ 



J2p-i{V2gt) = V2g5i^p. 



(4.30) 



Since the kernel of the reduced system f l4.23p is the same as the BES one (13. 3p . it is possible to 
use the same procedure introducing the functions 



+00 +00 

a^^\V2gt) = J2 Si';\9) J2p{V2gt), a^^\V2gt) = J2p-i{V2gt) 
p=i p=i 



(4.31) 



together with 



Sil\g) = 2{2p) I ^aX'{t)J2p{t), ^^;^i(^?) = 2(2p - 1 



dt 



Ik) 



(k) 



dt 



j(r'_\t)J2p-iit). (4.32) 



And, from the system (I4.23p . we derive the following equations for the functions a^\t): 



"+°° dt 

T 



af{V2gt) a^^\V2gt) 




J2p{V2gt) 



+00 ^2A;-1 

dt . ^ . ,^ J2p{V2gt) 



J2p-i{V2gt) = — 



sinht/2 



+00 ^2A:-1 



(4.33) 



dt- 



Att Jo sinht/2 



J2p-i{V2gt). 



The next step is to perform some manipulations on systems (14.30t I4.33p . in order to exploit their 



similarities. Concentrating ffist on (I4.30p . we multiply both sides of the ffist equation by Si^p'ig), 



®The use of a'f ^y^gt) is redundant with respect to S''^°^(fc) (|3.2p . However, since in Appendix A we will use 
results of [TB], we prefer to use here notations of [TO] . 
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and both sides of the second equation by S2p^_i{g). Summing over p in both of them, we end up 



with 







dt 
T 

dt 

T 




0, 

V2gS^^\g), 



where we notice that the coefficient Sf^\g) is exphcitly singled out. 

The same procedure can be repeated upon (14.331) . by multiplying the first equation by S^{g), 

AO) 

'a^^\V2gt) ai'\V2gt) crf{V2gt) a['\V2gt) 



the second by S^_i{g) and finally summing over p. The result is as follows: 



+°° dt 
T 

-"^di 
t 



1-e- 



+00 



— / dt 

in 



t 



sinht/2 



a 



(0) 



a'y{V2gt) a[!:\V2gtl ^ a^{V2gt) a'^\V2gt) 



1-e- 



e* - 1 



— / dt 

An 



t 



sinht/2 



a 



(0) 



iV2gt) 



iV2gt) 



A direct comparison with the previous equations allows to eventually obtain the integral represen- 
tation for S[''\g), 



V2gSPig) = -^ 



An 



+00 ^2fc-l 

dt- 







sinhi(;/2 



af\V2gt) - a^°\V2gt) 



(0), 



(4.34) 



For what concerns the coefficient Sl{g), relevant for the computation of fi{g), the procedure is 
identical - one starts from (13. 9p - but the result is slightly different. We have 



V2gS['\g) 



dt 



t 2cosht/4 



e--^(r^^\V2gt) + e^af\V2gt) 



(4.35) 



Equations (I4.34[ I4.35P are the representations of, respectively, S[''\g) and Sl{g) in terms of the 
BES density. 

Now, for completeness' sake and since we will be using it, we write also the integral representa- 
tion in terms of the solution of the BES equation for the density and its derivatives in zero a^^^'^'^\ 
Such representation trivially follows from definition (13. ip : 



-ASsfi + I'^dsa^H 



(0);(s) 



i ^dsCF^ 



(0);(s) 







+00 IS °° 

dk-^^J2^p^(9)Jp{V2gk) 



sinh k/2 



(4.36) 



+ 00 



dk 



p=i 
" .(0) 



sinh k/2 



(rf{V2gk) + (r^^\V2gk) 



We are naturally arrived at one of the main points of our work and some comments are now in 
order. We have found that the generalised scaling functions fn{g),n>2 enjoy an expression (14.221) 
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in terms of the 'reduced' coefficients 5*} (g) - related to the solution of the BES equation - and of 
the densities and their derivatives at zero a^'^''^'^^\ n' < n — 3 - related (when n' > 1) also to the 
other, 'higher', systems of equations (13.91) . (I4.23p . 

From the physical point of view the quantities S[''\g) and 

^(n');(s) g^^g 

expressed in terms of 

the various masses of the SYM theory. In next section, we will show that all these quantities can 
be explicitly computed in the strong coupling limit by using iterative methods. In such a limit 
they will be all expressed in terms of one mass, the mass gap of the 0(6) nonlinear sigma model, 
embedded in the Af = 4 SYM theory. 

5 Explicit results at strong coupling 

In this section we will compute the strong coupling limit of two sets of quantities: on the one hand 
cr{o);(s) g^^^ sl'^^g)^ on the other hand 0"^"'^'*^*^ with n' > 1, and c*^™\ 

The ffist set of quantities can be computed by relying on the solution of the BES problem. 
Indeed, a simple but long calculation, which was done using the method developed in [16] and 
whose details are given in Appendix A, allows to give the following analytic estimates at strong 
coupling, 

V2gS[%) = ^P@"2.n(,) + 0(e-^), 

ap-''-''^ = -{^y rrm{g) + 0{e~'^), k>0, (5.1) 
^(o);(o) ^ 4-7rm((7) + 0(e~'^), 

where m{g) (13.131) agrees at the leading order with the mass gap of the 0(6) NLSM. Using the 
ffist of (14.361) we can write, for the all loops density: 

V2gSi'\g) = t^(|)''2m(^) + 0(e-='^), 

(5.2) 

aC^W^'^) = -[^y\m{g) + 0{e-'^^). 

For what concerns the cr*^"'-''*-*'', with n' > 1, and the c*-™\ we ffist solve at strong coupling the 
system (I4.23P for the Sr''\g), finding an expression for such coefficients as asymptotic series in 
inverse powers of g. Then, using this expression, we will write a recursive equation relating the 
large g leading behaviour of the quantities a^"'-''*-*'', with n' > 2, and the constants c*-™-*. This 
equation - together with (I4.27p - will allow to find recursively both cr('^')'(''), with n' > 2, and the 

Having exposed the plan of our work, let us start from the reduced system (14.230 and let us 
look for a solution of it in the form 

n=k ^ n=k+l ^ 
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with 

S2r^ = 2m— -TT -IK , n>k, 5.4 

i (m — n + 1) 

S*iT" = (2,n. -!) [:;'"+" -'' (-ire., ">i + l. (5.5) 

L [m — n + 1) 

Usual techniques [151 [Hi allow to find the unknown ^2n-i)^2n solutions of the two recursive 
equations 

n—k r(fc) 
m=0 



(5.6) 

,_i\n 92n-2fc+l _ 1 n+l-fc r(A;) 

2"+^ 2n {2n-2k + 2)\ ^" ^'^+^ ^ ^ (2m)! ' 



m=0 



By comparing such equations with the corresponding equations for the coefficients appearing 
in the asymptotic solution for S]^\g) we find the simple correspondence 

b^N=-^ bN-2k, N>2k. (5.7) 

Putting all the relevant relations inside (14.22^ and redefining (for conciseness' sake) the indexes I 
and s, we finally find the asymptotic expansion^ 



k=l 1=1 s=0 ^ ^ 



^ 2'-(-l)"' T{p + n' + l-s) 1 ^ 'jl' {c^^)y- 

2^ 2n'+2Z-2. r(p-n'-/ + S + l) '"'J ^ J-J- J^! ' ^ 



and 



I — 1 ;— 1 o— n \ / 



s{n-k)-l2s) 



k=\ l=\ s=0 

2'-^(-l)"' r(p + n' + Z-s), 1 Vf'(c^^ 



where the positive integers jm have to satisfy 

A;-2Z k-2l 



J2jm = 2l + 1, Y^mj^ = k. (5.10) 



m=l m=l 



^The notation [x] present in (|5.8l 15. 9p stands for the integer part of the semi-integer x. 
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We are now ready to discuss the strong coupling behaviour of the densities of Bethe roots and 
their derivatives at m = 0, 0"*-"'''*^'*^ when n > 1. We have to distinguish between the case n = 1 
and the cases n > 2. 

In the case n = 1 we have that 



0", 



{i);(o) 



-4 In 2, i-'dsa^o^'^"^ 



(4-2^+')r(s + l)C(s+l), s>2, 



for the one loop theory and 



(!);{«) 

H 



~00 7.S 



sinh k/2 ^ 



(5.11) 



(5.12) 



for the higher than one loop contributions. We now insert in (15. 12^ the asymptotic solution [15] for 
Sp^\g), reported also in (13.101 [3?TT]) . Summing first over p and then over n' and finally integrating 
in k we get 



a 



(i);(o) 

H 



TT 



(i);(^) 



31„2--, 

(2^+2 _ 4 + 2~2^ - 2~')T{s + l)C(s + 1) 



TT 



s+1 



s>2, (5.13) 



where Ek are the Euler's numbers. Summing with the one loop results, one gets the explicit formula 



a 



(i);{0) 



-ln2 - -, 
2 ' 



?-^4cr(')'W = (2-2« _ 2-')r{s + l)C(s + 1) - 
For the case n > 2 we start from (13.141) and write 



1^. 



s > 2. 



dk . 
— I 
27r 



Tck 
sinh 



-S^''\k) 



e 2 
sinh 



dp ,sin(A; — j9)c 
or(p) 



27r 



k — p 



(5.14) 



(5.15) 



2 """" 2 

We then insert (15.81 15. 9p in the expansion (13.171) of S^'"'\k) in series of Bessel functions and use 
(14.201) to express the integral term in the square brackets. Summing on the indexes p and n' coming 
from (15.81 15.91) . we end up with the previously announced recursive equation, 



n [ 2 J / o; \ 

E E M)T ( 2!' ) • 

k=l 1=1 s'=0 ^ ^ 



+~ dkk'+^'-^'' 



71 sinh I \ cosh k 



k 

62 _fe 

— e 2 



fc--2/ 



E n 

{ii,-jfe-2j"^=i 



jrn- 



n [ 2 J / 0/ \ 

E E M)T ( 2,' ) ■ 

fc=l «=1 s'=0 ^ ^ 



(5.16) 



^2-s-2i+2s' _ 2-2s-Ai+is'^Y{s + 21- 2s' + l)C(s + 21- 2s' + 1) + 



+ 



T^x s+2/-2s'+l 



i) 



fc-2/ 



+2i-2s' 



E n 

{ii,-jfc-2i}™=i 



Jm- 
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where, again, 

k-2l k-2l 

^jm = 2/ + l, '^mjm = k. (5.17) 

m=l m=l 

As we said before, such equation has to be solved together with f l4.27p . 

To summarize the results, in this section we have shown that, similarly to the case of Sr^\g) 
|15j . also the system for the Sr''\g) can be solved at large g. This allowed to write for the cr'^"')'*^'^^ 
a further recursion relation f l5.16p - which goes together with the explicit expressions (15.141 15. 2p . 
coming from the solution of the systems for si^\g) and from results of [12], respectively. In this 
regime, the set of constants c*^™'^ and the set of densities o"^"'^'^*^^ can be computed by solving 
simultaneously the recursive relations (I4.27P and (I5.16p . Putting their expressions, together with 
the solution for S[''\g), into Km . one will get the expression for s["'\g) and, consequently, for 
fn{g) = V2gS'C\g) at strong coupling. 

As an application of all these techniques, in the next section we will compute explicitly the 
strong coupling limit of the scaling functions f^^g), . . . , fs{g)- 

5.1 Examples: fs{g) to fs{g) 

The previous machinery can be tested by computing the strong coupling behaviour of fn{g), for 
3 < n < 8, in order to compare it with the available result from the 0(6) NLSM. 

First of all, we need to know the expression of c^^\ . . . , c^^\ From the recursion formula (14.270 

we get 

'''' = (5-^8) 

^ - y [(j{0);(0)]4 ^[^(0);{0)]3 ' ^^■'^'^> 
^(3);(0) 2 y0);(2)^m0) r^(l);(0)l3 ^ ^(1);{2) 

" [^(0);(0)]2 3" [^(0);{0)]5 ^ " [^(0);{0)]4 ^ g " [^(0);{0)]4 ' ^^"^^^ 

_ TTfqW'W]^ 5^3^(o);(2)[^(i);(o)]2 ^(i)'^^) ^(1)' ^ 27ra(3);(o)cx(i);(o) 

^ ~ [o-{0);(0)]5 3[cr(0);(0)]6 3[o-(o);(o)]5 \^iomf ^ 

^5[^(0);(2)]2 TT^gW-W 7raW'(o) 

+ 1 onr-rnvrmifi + TrTrvTTnTW ' (5.22) 



12[a(o);(o)]7 120[a(o);(o)]6 [a(o);(o)]2 



(6) _ vr[a(i)'W]5 107rVW'(2)[a(i)'(o)]3 57r3a(i);(2) [a(i);(o)]2 37ra(3);(o) [a^i)'^]^ 



[o-(o);(o)]6 3[cr(o);(o)]7 3[o-(o);(o)]6 [o-(o);{o)14 

7^5[^(0);{2)j2^(l);(0) ^5^(0);(4) ^(1);(0) 2710^^^'^^^ a'^'^'^'^'^'' 

12[a(o);(o)]8 20[a(o);(o)]7 ^(oMp 6[a(o);(o)]7 



+ 



120[aW;(0)]6 3[(t(0);(0)]5 5[^(0);(0)]4 [^(0);(0)]2 



(5.23) 
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Then, referring for the notations to (14.201) . we have 



+ 



2 / 


■°° * ir fh 


— 






* lUH 

-00 2vr 




- r4(p,p)] 


2/ 


-00 2vr ^'^^^ 


-p,p) 


- r5(p,p)] 


^/ 




-p,p) 




4 7r" 


3[^(l);{0)]3/^2 


1 tt' 




3 


[o-(o);(o)]5 


15 [ 


cr(0);(0)]5 


^/ 


-00 2vr 


-p,p) 


- r7(p,p)] 




j^(l);{0)]4/^2 


^7[^{0);(2)]2;^2 



2 ,/iV(i)'W 



3 [cr(o);(o)]3 

TT^/l^ ^3[^(l);{0)]2/^2 ^5^{0);(2) /^2 



(5.25) 

(5.26) 



60[(t(0)'(0)]4 [o-(0);(0)]4 15[^(0);(0)]5 

^5/^4^(l);(0) 2 7rV(3);(0)/,2 ^ ^5^(1);(0) ^(0);(2) /^2 

15[a(o);(o)]5 ~ 3 [a(o);(o)]3 + 3 [ct(o);(o)]6 + 

(5.27) 



4 



2520[(t(o)=(o)]6 6[a(o)'(o)]6 ^ 126[(t(o);(o)]7 + 



,2 



3[(t(o)'(o)]6 ^ 36[(t(o)'(o)]8 [(t(o);(o)]7 420[a(o);(o)]7 ^ 3[a(o)'(o)]6 ^ 

2^3^^(l);(0)]^(3);(0)/^2 27rVW'W/l2 

+ [ct(o);{o)]4 3[a(o);(o)]3 ^^-2^) 

r dp , , , 7r7[a(l)'(0)]/i6 ^5r^(l);(0)l3/^4 ^7^(0);(2)r (l);(0)l/^4 



00 



27r ' "^"^^ ^'"^"^'^ 420[a(0);{0)]7 3[a(o);(o)]7 18[(t(o)^(o)]8 



7rV(^)-(2);,4 7rV(3);(0)/,4 27r3[a(i);(o)]5/^2 y^5^(o);(2) [^(i);(o)]3/^2 27r^[aW-(2)]2 [a(i);(o)]/^2 



TT 



7^{0);(4)[^(l);(0)]/^2 ^5 [^(1);(0)] 2^(1);(2) /^2 ^7^(0);(2) ^(1);{2) /^2 ^7[^{l);(4)]/^2 



60[(t(o);(o)]8 [a(o);{o)]7 18[a(o);(o)]8 420[a(o);(o)]7 

4^3[^(l);{0)]2^(3);{0)/^2 ^5^{0);(2) ^(3);{0) /,2 ^5 [^(3);{2)] /^2 27r3 [(T^^); 

[cr(0);{0)]5 ^ 3[^(0);{0)]6 15[cr(0);{0)]5 [(j(0);(0)]4 
2^3[^(5);(0)]/^2 

3[(t(o);(o)]3 ■ 



(5.29) 
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This implies that for /^{g), . . . , fs{g) we can give the exact (i.e. vahd W g) expressions: 



fsjg) 

2\p2g 



-TT 



6 [a(o);(o)]2 



h{9) 
2V2g 

hjg) 

2V2g 

hjg) ^ 

2V2g 30[a(o);(o)]^ 



1 a^m^'^ 

^ Z_ Q^^l ( \ 

[(X(0);(0)]3'^l ^3), 



120[a(o)'(o)]4 



2 

1 7rV(3);(o) 
'3 



30[a(o)'(o)]5 



+ 



(5.30) 
(5.31) 
(5.32) 

(5.33) 



+ 77 



1 ^5^(l);(0)^(0);(2) ^ 3 TrV^^^ ^ j ttV^K^) ' 
6 [(X{0);(0)]6 3 [^{0);(0)]5 " ^ [ct(0);{0)]5 



2v^(7 5040[a(o);(o)]6 



TT 



[^(i);(o)]^ 



+ 



^7^(0);(2) 



12[a(o);(o)]6 252[(t(o)'(o)]7 



(5.34) 



+ 



+ 



57r3[a(l)'W]^ ^7[^(0);(2)]2 ^5^(0);(2)[^(l);(0)]; 

6[(t(o);(o)]6 ^ 72[a(o);(o)]8 



2[(t(o);(o)]' 



^5[^(l);(0)j^(l);(2) ^3 [^(1);(0) j^(3);(0) yrV^^^'^") 



6[a(o);(o)]6 



840[(t(o);(o)]' 



+ 



+ 



840[(t(o)'(o)]7 
7r5a(3);{0) 



Cr(0);(0)]4 3[^(0);(0)]3 
^5[^(1);(0)]3 ^ ^7^(0);(2)^^(1);(0)] 



7^(1);(2) 



TT cr 



252[a(o)'(o)]' 



+ 



(5.35) 



3o[a(o);(o)]E 



^3[^(1);(0)]5 ^ 77rV(0);(2)[a(l);(0)]3 ^ ^7j^(0);(2)]2^^(l);(0)] ^ 



cr 



(0);{0)]^ 



9[^(0);(0)]£ 



+ 



^7^(0);(4)[^(1);(0)] ^5 [^(1);(0)]2^(1);(2) ^7^(0);(2)^(1);(2) 271^ [a^^^'^O)] V^^);^ 

120[a(0)^(0)]8 2[aW;(o)]7 + 36[aW;(o)]8 ^MiWp 

^5^{0);(2)^(3);{0) ^5^^(3);{2)j ^3^^(l);(0)j^(4);(0) ^7^^{l);(4)j ^3^^(5);{0)j 



30[(t(0)'(0)]^ 



[o-(o);(o)]4 



M0[(tW'(0)]7 3[^(0);(0)]3 



As stated before, there formulae express the /n(5')'s in terms of the masses of the TV = 4 SYM 
theory. 

In the limit g oo, the quantity cr^"-''^*^ for n > 1 can be computed using (15.140 and solv- 
ing together the recursive equations fl5.16p and fl4.27p . In particular we have, for their leading 
contributions. 



r(i);{o) 



TT 

ln2--, a 



(i);{2) 



3C(3) + TT^ 



(3);{0) 



TT 



(i);{2) 



^(i);(4) _ _Afvr5 I q/^fc:^^ ,(3);(2) _ 



(i);(4) 



cr 



(4);{0) 



6[(t(o);(o)]2^ 

^^(1);{0)^(1);(2) 



-TT 



3[ct(o);(o)]= 



(5.36) 
(5.37) 
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a 



(5);(o) 



^2^(l);(2)r^(l);(0)l2 



^4^(0);(2)^(1);(2) 



;(4)] 



2[a(o);(o)]4 30[a(o)'(o)]5 ^ 120[a(o)=(o)]4 ' ^^'^^^ 

These data, together with the hmiting leading values at strong coupling of a*^"^'*^*'' and S'^\g) fl5.2p . 
allow to obtain for fz{g), • • • , f%{g) the following leading values as (yf ^ oo, 



/3(^7) 



TT 



2Am{g) 



I2[m{g)]^ 



5l 



TT 



640[m(^)]3 8[m(^)]3 



53 

90 



120 



+ 



7182[m(^)]5 192[m(^)]5 2A[m{g)]^ I8[m{g)]^ 



+ 



+ 



m[m{g)]^ 16[m(^)]6 A[m{g)]^ 560[m(^)]6 6[m(c/)]6 280[m(^)]6 
where we used the compact notations: 

1 1 



(5.39) 
(5.40) 
(5.41) 
(5.42) 
(5.43) 
(^.44) 



2s+l 



TT 



— Ti-r 



n=0 



2s+l 



+ 



n + 1 



\2s+l 



(5.45) 



For instance, we have 
5i 



71 2 



5, 



47]-3 



[3C(3) + 7r3], ^5 



487r5 



[9C(5) 



TT 



(5.46) 



After a lengthy but straightforward calculation it is possible to show that such expressions agree 
(at the leading order) with the corresponding formulae computed in the framework of the 0(6) 
NLSM, i.e. the coefficients 2"'~^Qn given by the general formulae of 



5.2 Numerical evaluation of the next to leading correction to the mass 

In order to check results presented in the previous sections we numerically estimated (via a best 
fit procedure) the leading non-analytic contribution at strong coupling for various quantities: 

"i = /i(fl') + l, "2 = , as = -4 (5.47) 

a, = -[V2gS['^], a, = --^[V2gS?^]. (5.48) 

All of them turned out to be proportional to the leading 0(6) mass and the actual value was in 
perfect agreement with the computation by Alday and Maldacena. As a result, we were naturally 
led to conjecture the exact value for the leading exponential behaviour. 
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1 


2 


3 


4 


5 


1 





-0.818(2) 


-5.72(3) 


1.64(1) 


16.3(2) 


2 


0.818(2) 





-4.96(1) 


2.46(2) 


17.2(2) 


3 


5.72(3) 


4.96(1) 





7.35(1) 


22.04(15) 


4 


-1.64(1) 


-2.46(2) 


-7.35(1) 





14.7(2) 


5 


-16.3(2) 


-17.2(2) 


-22.04(15) 


-14.7(2) 






Table 1: Values of the differences 6i^j. The index i runs along the rows, and the index j runs along 
the columns. 



Furthermore, we are in the position to give a quite precise estimate of the next to leading 
correction to the mass m{g), which is known to be of the form 

k n 95/8„l/4 

m{g) = kg'/\l + J + --- k = (5.49) 

The analysis of the data at our disposal gives access to the quantity ki, whose best fit estimate 

is 

fci = -0.0164 ±0.0005. (5.50) 

It is interesting to point out that such a correction has to be the same for all the functions 
a[s at leading order in the strong coupling limit. We checked that it is actually the case, giving a 
further proof of the reliability of our computations. 

As matter of facts, all the corrections in 1/g, the Brezin-Zinn Justin function [22], should be 
the same for the previous observables, because the computations of appendix A show that the first 

- — a 

discrepancies should appear at order 0(e ). 

Such a fact can be checked numerically, by performing a numerical analysis of the differences 
Aij = ai — aj, i < j. With a best fit procedure, we have been able to verify that all the A's 
actually behaves as [fn{g)]^/g but with different amplitudes, namely 

Ay = *.,fcM! + ... (5.51) 

A consistency check that is fulfilled the numerical amplitudes 6ij is that the relations which are 
not linearly independent are actually compatible within the numerical precision. The amplitudes 
are collected in table [TJ 

6 Summary and outlook 

In this paper we have discussed a systematic method which aims at the computation of all the 
generalised scaling functions fn{g) appearing in the sl{2) sector of A/" = 4 SYM in the scaling limit 
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In our approach the generahsed scahng functions are expressed, for all values of (7, by fl4.22p . 
All the ingredients entering such expression can be found on the one hand from the solution of the 
BES equation and on the other hand from the solutions of the systems (13. 9 p for Si^^ (g) and (14.230 
for the reduced coefficients Si-''\g). At strong coupling the explicit solutions to the systems for 
the reduced coefficients has been explicitly worked out (15.81 15. 9p at the leading order. Together 
with previous analogous results concerning the BES equation [12j and Sl^\g) [T5], this eventually 
allowed to compute the strong coupling behaviour of the first eight generalised scaling functions 
(15. 39115. 441) . Our results agree with the analogous findings [H] coming from the 0(6) nonlinear 
sigma model. 

For what concerns future work, one has to say that in the s/(2) sector of Ai^ = 4 SYM other 
regimes - e.g. the large j regime in the context of limit (11.41) [23] - and, more generally, limits 
different from (11.41) are relevant for comparisons with pure string theory results. In this respect, 
the limit s, L — >■ 00, g 00, / = L/{g\ns) fixed - the so-called 'semiclassical scaling limit' - has 
been widely studied [24j. Application of our techniques to this case is a possible future direction 
of investigations. 

Finally, one has to mention the new line of research related to the recently discovered duality 
between J\f = 6 super Chern-Simons (SCS) theory with U (N) x U{N) gauge group at level k and 
superstring theory in the AdS4 x CP^ background, when is large and the 't Hooft coupling 
A = N/ k is kept fixed [25] . Integrability on the gauge side [26l [27] and on the string side of the 
duality [28] was shown. Bethe ansatz-like equations were proposed [29] for the SCS theory and 
tested in various ways [30] . It could be surely of interest to apply the techniques discussed in this 
paper also to this new field of activity. 

Acknowledgements We thank F. Buccheri, D. Bombardelli, F. Ravanini for discussions and 
suggestions. We acknowledge the INFN grant "Iniziativa specifica PI14", the international agree- 
ment INFN-MEC-2008 and the Italian 'PRIN' for travel financial support. The work of P.G. is 
partially supported by MEC-FEDER (grant FPA2005-00188), by the Spanish Consolider-Ingenio 
2010 Programme CPAN (CSD2007-00042) and by Xunta de Galicia (Conselleria de Educacion and 
grant PGIDIT06PXIB296182PR). 

A Non-analytic terms at strong coupling 

The aim of this subsection is the analytic computation at large g of the leading contributions to 
eqs. (I4.34P by means of the techniques developed in [T6] . 

In particular, we will calculate the large g behaviour of the following integrals. 
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First of all, we make use of the BKK transformation [121 fT6] . 



2af(t) 

in order to rewrite the integrals as 
B2k-iig) = 

C2k{.g) = 



dt ( t 



t \V2g 



2k 



Sf (t) itanh^E^^^m 
\/2 ^ ' ' 



cosh 



(A.l) 



sinli — - — Lvwoxi 

- (t)) + (t)) 

,cosh^ cosh^ 



+00 



dt ( t 



t \V2g 



2k 



cosh -7=- sinh 



The BES equation can be rewritten in terms of the functions T,^^^ [16], with |m| < 1 

r+00 

/ dt sm{ut) (t) + (t)] = , 
Jo 



2V2g ,^{0) 



(A.2) 



+00 



rft cos(Mt)[sL°^(t) - Sf (t)] = 2(2^2(7) 



and the ratios of hyperbolic functions admit a useful integral representation 



sinh 



t 



2k-l 2V2g 

sinh 



+00 (i^A;-! 

COs(-ut) — ri — 7 



2k 2V2g 



(-1)'' g / ciu sin(Mt) 

J —00 

r+00 

{-!)'' g / (i-u sin(M)f:) 



du'^'^ ^ 



cosh \plg'KU 



cosh^ 



+ CXD 



COs(Mt) 



dxi^^ 

d^ 

du^^ 



cosh v^(yf 
cosh \/2g 



TTU 



TCU 



Plugging them in the integrals -62^-1, C2k we obtain 

2fc /.+00 



B2k-i{g) 



'-^fg 



1 



+ 



sin(ut) 



—00 



(it cos(ut) 



2fc-l 



cosh V^fyfTTM 



sinh^ 



1 



cosh \plg'KU 



C2k{g) 



-i)'g 



V2g 



2fc+l ^+00 



du 



(sL°^(t) + sf (t)) 



cosh A/25f 



(0), 



(it cos(nt) 



(iw^'^ 



1 q-KU 

cosh^ 
cosh ^J^gnu 



(sL°)(t)-sV'^(t)) + 



(0), 



+ 



+00 ^2fc 

dt sin (-ut) , „■ 



cosh \plgTXU 



(sL°^(t) + sV"^(t)) 



(0), 
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Let us evaluate them in the large g limit. The strategy is to split the integral over u in two 
intervals |m| < 1 and 1^1 > 1 in order to use the constrains (1A.2I) . The former gives, together with 
the use of such constraints: 



1 c;2fc-i 
du 

1 

1 



2{2V2g) / du 



cosh \/2gTxu 



+ 00 



dt cos(Mt)(sL°^(t) -SV'^(t)) 



(0), 



du 



1 
1 

-1 



du 



du^'^"^ 

d^ 
du^^ 



du^^ 



cosh \plgTxu 



cosh \plgi{u 



cosh^ 



+ 00 



rft sin(Mt) (EL°^(t) + EV"^(t)) = 



(0), 



cosh ^/2g7ru 



+ 00 



dt cos(Mt)(sL°^(t) -Sf (t)) 



(0), 



I du^^ 



cosh \/2g'KU 



^ d^^ 
du ■ 

1 



du^^ 



cosh \/2gTiu 



dt sin(Mt) (sL°^(t) + SV'^(t)) = 0. 



(0), 



Since we are interested in the large g behaviour, we can perform the previous integrals by rewriting 
them as the difference of the integrals with support over (— oo, +cxd) and (— oo,— 1), (1,+cxd), and 
finally taking the leading exponential in the integrands, so we will have for n > (we will use a 
single index n because at this order there is no distinction between even and odd indexes): 



A{2V2g) 



du - — e 



8V2g 



Tig_ 
V2 



n-l 



e"^ + 0(e"^) . 



The case with n = needs to be treated separately, because we also have to take into account the 
contribution of the integral over (— cx), +oo)\ 



du 



cosh^ 



cosh \/2g 



TVU 



+ 00 



du 



cosh ^j2gTXU 



+00 



du 



cosh \/2gnu 



We end up with 

2{2V2g) 



du 



-1 



cosh \f2g'nu 



aV2 - 8V2g I ^ 1 e"^ + 0(e"'^) 



(A.3) 



We stress that the above integrals have the same structure for any n only at leading order for large 
g. Next to leading orders will differ because of the different form of the integrands. 

It is possible to show that, taking the leading exponential only, the integrals over l^l > 1 take 
the form 



g-KU 

due~^ 



/ cos(Mt)(sL°^(t)-Sf (t))+ / dtsin(Mt)(sL°Ht) + sf (t)) 

Jo Jo 
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which is accurate up to 0(e ^) terms. The integral was estimated in [16], and taking into 



account the difference with the notations of that papecj, we have 



due 



+00 



dt cos(ut) (EL°^(t) - (t)) + / dt sm{ut) + (t)) 



L^o 



(0), 



TT 



—= m(q) H — e ^ 
V2 vr 



(A.5) 



where we pose, as before. 



m 



{g) = kg'/'e-^^'[l + 0{l/g)], k 



25/8^1/4 



(A.6) 



It is interesting to notice that the 0(1/ g) in the previous equation stands for power-hke corrections 
which in principle can be computed by inspecting the sub-leading terms of the Ihs of ( lA.Sp . 
If we put everything together we have, up to the level of accuracy discussed before 



B2k-i{g) 
C2k{g) 
C,{g) 



16v^ /7r\2fc 2V2 / 



TX 



2k 



TT \2 

n\ 2fc 



TT 



— m{g) + -e ^ 

V2 TT 



8V2 _^ 



TT 



_9ZL , /- ( TT , ^ 8 _92l\ 

e v^+4 + V2 -= m(q) + -e ^ ] . 

V V2 ^ / 



— m(^) + -e 
V2 TT 



We notice that the first and the last terms always cancel, hence we have eventually: 



^2k~i{g) 
C2k{.g) 



TT 



2k 



\k+l 



TT 



2m{g) + 0{e Vi) , 
Tim{g) + 0(e~ 



> 



C,{g) = 4-7rm((7) + 0(e~'^), = 



J0);(0) 



Therefore, we end up with the following estimates at large g (i.e., up to the order 0(e v^)): 



^9~sf\g) = -^B2k-i{g) 



\k+l 



An 



a 



{0);(2fc) 
H 



in 

2k 



TT 



2k 



2) 2m((7) 



a 



(0);{0) 
H 



i C2k{g) = - [j^) nm{g), k> 
Co{g) = A - TT m{g) . 



(A.7) 



^"it is easy to work out the following relations between our g, and the same quantities gsK, in the paper 



9 T.(o) ^± 



(A.4) 
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We conclude with few words on the coefficient relevant for the computation of fi{g). Starting from 
(14.351) and using the BKK trasformation (lA.ip . we get 



+00 







dt 
2i 



sinh 



- cosh 



^(sL°)(t)-sf(t)) + 

V2g 



+ 



(A.8) 



which, after the use of the map flA.4p . coincides with (53) of [16j. Therefore, the strong coupling 
analysis can be done following that paper. 
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